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A New Approach to the Riemann Hypothesis 


І. C. Chakravartty 


ә 


1 Introduction 


The Zeta-Function of Riemann is defined by the Dirichlet series: 


C(s) = Уп, 
n=] 


where s = ø + it, с and t being real numbers. 

The series is convergent and the function is analytic for o > 1. It 
is a meromorphic function having a simple pole at s = :1 with residue 1. 
Although this function was known to Euler it was Riemann who discovered 
its most remarkable properties. Since then, this function has proved itself 
to be of fundamental importance, not only in the theory of prime numbers, 
but also in the higher theory of Gamma functions and allied functions. 

It was first conjectured by Riemann [2] in 1859 that the complex (non- 
trivial) zeros of the Riemann zeta-function in the “critical strip” 0 < ø < 1 
lie on the “critical line” o = 1/2. This conjecture has now become famous 
as the “Riemann Hypothesis”. It is highly probable that Riemann's con- 
jecture about the non-trivial zeros of the Riemann zeta-function is correct, 
but it has not been proved or disproved to this day. This famous hypothesis 
is now more than one hundred and thirty-five years old and many world- 


renowned mathematicians of successiye generations devoted a considerable 


/ 
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&mount of time and energy to resolve this hypothesis, without &Dy Success 
so far. Some prominent mathematicians working on this subject now be- 
lieve that there has to be some new approach to this problem. The object 
of this paper is to suggest a new approach to this age-old problem which, 
I hope, will someday lead to a solution to the problem. 


2 Some Useful New Notations 


It has long been known that there is a close connection between the prime 
numbers and the non-trivial zeros of the Riemann zeta-function. Riemann 
himself noted the existence of such a formula [5]. One such “explicit for- 


mula", given by Von-Mongoldt, is as follows [6]: 





$(z)—z = Y A(n-z 
2 узае аа C), 
= исе 2n EO)" (=> 1) (1) 


In this formula, the notation У is used to indicate that the term n = x 
if it arises, is to be halved. Also, if n is a positive power of a prime, say 
п = р", p being a prime and m bene а positive integer, then: A, = log p. 
Otherwise, A, = 0. Finally, p runs through the complex zeros of ¢(s) and 
dp із the principal value of: 
pm, 
—T«Im(p)«T 

Thus, the formula (1) connects & sum involving powers of primes and the 
imaginary parts of the zeros of the Riemann zeta-function. 

In almost all standard texts and papers written on the subject of the 
Riemann hypothesis, the symbol p is used to denote the complex zeros of 


the Riemann zeta-function. This convention, however, has proved to be 
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a disadvantage when it comes to getting a good insight into the problem 
of solving the Riemann hypothesis. Furthermore, I do believe that this 
convention actually hindered research in the subject. Instead of using p, 
we сап, to our great advantage, use 1/2+iy to denote a complex zero of the 
Riemann zeta-function, irrespective of whether the Riemann hypothesis is 
true or false. If the Riemann hypothesis-is true then all the values of y in 
1/2-F iy are real numbers. On the other hand, if the Riemann hypothesis is 
false, then some of the values of y in 1/2+ iy will be complex numbers. In 
either case, since complex zeros of ¢(s) occur in conjugate pairs, the series 
— x = On the right-hand side of equation (1) becomes: 
pitty gi 
Еа = +e) (2) 
S i + у 5 — Vy 
lo make this simpler from the viewpoint of calculations, we can replace 
1/2 + ty by +iy and z by e7. The expression in (2) then becomes: 





туту —t7 
XLV 
ji d 
Е y sin yy 








Se mk CH yee „у, @ 


p P7k f£ Qn+k C(k)’ 
Putting k = 1/2 and a = еу, the first term on the right-hand side of (3) 
becomes: 


$ Р`?“1орр= У pi logy. 
p" «ey mlogp<y 
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The second term on the right-hand side of (3) becomes: 


= 7 (anri) 


2п +$ 


1 х 
1-2 zi 


907% + tan! 27 — 5 log 


n=. 


= —2e1 + tan”! e às — z log tanh y 
‚ We now consider the functional equation of the Riemann zeta-function: 
1 
C(s) = 281°! sin 557 T — s)C(1 — s). 


Taking the logarithmic derivative of this equation and then putting s = 
1/2, we obtain: 

С) q. 4. 3 

— L = С — -m- r-log87, 

c) 2 "ME АШ 
where C represents the Euler's constant. We thus obtain the following 


result from (3): 


Уу] pt logp – 2e? 


mlog psy 
sinyy 1 1 1 
= 932 — -C – -r — –1ор8т – 2 
b E 2 47 — 5 log8r e 
1 
-iy + tante- — 2 logtanh cy, (y > 0). (4) 


3 The Prime Number Summation Formula 


The question of finding a “general”, rather than an “explicit formula” , con- 
necting the non-trivial zeros of the Riemann zeta-function and the prime 
numbers have been the subject of research for a long time now. It was 
pointed out by Winter [8| in 1943 that this relationship appears to be of 
a reciprocal nature, analogous in some way to the “Fourier Inversion For- ' 
mula". In 1948, A. P. Guinand came out for the first time with a general 
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formula connecting the non-trivial zeros of the Riemann zeta-function and 
the prime numbers. This formula [3] actually connects the imaginary parts 
of the non-trivial zeros of the Riemann zeta-function and the logarithms 
of powers of primes as arguments of the functions concerned. With appro- 
priate conditions on the functions f(x) and g(x), this general formula can 
be stated as follows: 


| X jml- | five 


иша (LESE p 
1 f% 1 ей 
“af, f Е - e) © 
1 { 
= -У?т Y o) +- [ ot) log dt 


where, p runs through the prime numbers, m through the positive integers 
and p = 1/2+iy through the non-trivial zeros of the Riemann zeta-function. 
It must be noted in this connection that if multiple zeros of ¢(s) exist, then 
the corresponding terms g(y) in this formula are to be multiplied by the 
order of the zeros. 

Later, A. P. Guinand [4] discussed more general relationships of this 
type and showed that the relationship can be regarded as a type of Fourier 
reciprocity for almost periodic functions. In 1952, Weil [7] published a 
result for a wider class of zeta-functions. 

However, all these results were proved assuming the Riemann hypoth- 
esis to be true. The question naturally arises, can we establish a general 
formula connecting the zeros of the Riemann zeta-function and the prime 
numbers without assuming the Riemann hypothesis? In 1970, I was able 
to prove a general summation formula [1], which, unlike Guinand’s for- 
mula, was not based on the assumption that the Riemann hypothesis is 


true. This formula, therefore, is not dependent on any as yet unproven 
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hypothesis. With appropriate conditions on the functions ф(х) and x(x), 


this summation formula can be stated as follows: 








p of e i" 
-3], E — iv) o(u)du 
= -Va| Z v- ve edu), © 
27 Jo 2T 
where, ф(х) = v2 Jo Y(t) cos atdt, p runs through the prime numbers, m 
through the positive integers and 1/2 + iy (where some of the values of 


ү are complex since we aie not assuming the Riemann hypothesis) runs 


through the complex zeros of the Riemann zeta-function. 


4 A New Perspective on the Riemann 


Hypothesis 


We will now use a particular case of the Prime Number Summation Formula 


to gain & rather new perspective on the problem of solving the Riemann 


hypothesis. 
Let 
P) = T (e Mee 4 ehet) cos zy 
g(x) = e em соз 2(2 — y) + e dot T COS (т + y) 


With this pair, we find that the right-hand side of [5] can be split into: 


— Vin У, e M79" /T* cos 2(y — y) 
and —V2- Y, e- Mor /T^ cos z(y + y) 
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We, therefore, obtain the following result from the summation formula: 


M -{т -}(z-mlogp)?/T? 
Jim. 2. log p (e 


+e7a(e+mlogp)?/ E cos(y m log p) 


| т, іу = +5 


(7) 
0, elsewhere. 
Reversing the roles of f(x) and g(x), we obtain the following reciprocal 


formula: 


Т, —H(y-z)?/T? | 2 4(7+а)2/7° 
jim БРЗ (е pe ) cos(y7) 


_ | rlogp, ify = mlogp (8) 
0, otherwise. 


From (7) and (8), we see that there exists a Fourier reciprocity between 
the weighted sequences (m log p, pi" logp). 

Formula (8) remains true even without the Riemann hypothesis. How- 
ever, if the Riemann hypothesis is false, then the formula (7) needs modi- 
fication if 1/2 + iy is in the neighbourhood of an exceptional zero of С (s). 

If the Riemann hypothesis is true, then the weighted sequence {m log p, 
p i" log p) has just those fluctuations of spacings and wéightings (and only 
those) permitted by a sequence of real frequencies {7}. Consequently, any 
assumption that the Riemann hypothesis is false (i.e. some y’s are complex) 
should imply greater fluctuations in the former sequence. 

The question then naturally arises: Are there limits to these fluctua- 
tions? There are. The values of mlogp near z = logn can be spaced no 
Jess than the logarithms of consecutive integers near z. But if n is large, 


log(n + 1) — logn ~ 1 = e^*. Therefore, the minimal spacings near z is 


about e^. It seems likely that the Riemann hypothesis depends on the 
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limits of these fluctuations and that an investigation of these limits may 


some day lead to a solution of the Riemann hypothesis. 
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Some Results of M ^/M/1 Queue with 
Generalized Vacation 


Gautam Choudhury 


Abstract: The Laplace Stieltjes Transform (LST) for the joint 
distribution of the queue size and vacation period has; been derived 
for an M* /M/1 queue with vacation. 


Keywords: M* /M/1 Queue, Queue size and Vacation period. 


1 Introduction 


Recently, Borthakur and Choudhury (1997) have considered an Мх /G/1 
queue with generalized vacation, where the server is turned off each time as 
soon as the system becomes empty and he goes on vacation for a random 
length due to some uninterruptible task (e.g. machine breakdown, machine 
repair, scanning for new work etc.), which we shall refer to be a generalized 
vacation period. A generalized vacation starts at the instant the system 
becomes empty and the server leaves for the first vacation and then for a 
second vacation and so on until the server finds one or a batch of customers 
in the waiting line after returning from each vacation. 

Due to the theoretical structures and applications in number of real 


life situations such as inventory/ production systems [e.g. see Lee et. al. 


?AMS Mathematics Subject Classification, Primary : 60K25; Secondary : 05A15. 
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(1994) and Choudhury (1997)], digital communication and data transmis- 
sion systems [e.g. see Medhi (1997)], vacation queues become the subject 
of deep study in recent years. Hence, in this article we derive the LST of 
the joint distribution of the queue size and vacation period due to vacation 


for an M*/M/1 queue with generalized vacation. 


2 Main Result 
First of all let us define following notations: 
A - arrival rate, 
ш — service rate, 
X – arrival size random variable, 
be - R(X =k) 
bU) — j-fold convolution of by with itself, and bf) = 1, 
X(z) — Probability generating function (PGF) of X, 
V — vaction time random variable 
V(x) — Probability distribution function (PDF) of V, 
V*(0) — LST of V and 


p = (AE(X)/u)(« 1) is the utilization factor or offered of the system. 


Borthakur and Choudhury (1997) have obtained the joint distribution 
of the queue size and vacation period P,(x) (= > 0,j > 0) as а product of 


two independent random variables, which can be expressed as 


Рх) = L(z)A;(z); 22 0, j = 0,1,2,... (1) 


Some results of Мх /M/1 queue ... 11 

where I(x) = = ad “4 ів the PDF of the length of a vacation period 
and | 

2 ,—-Acr k 

A;(x) = Р, | arrivals during (0, z]] = У? e Аа) Ы, 

k=0 

Let P*(8;z) be the PGF of (P?(0);j > 0} [the LST of P;(x) ; x > 

0, j 2 0], then proceeding exactly in similar manner with Choudhury 

(1996) [see section 6], we get from equation (1) 


P*(6;z) = У`гіР"() (:|21<1:) 
7=0 


_ (1=p)fl-V*(0-+A~AX(2))} T 
E(V)[8 + A — AX (zj] 


Now if we take limit 0 — 0 in equation (2), we get 


lim P*(;2) = (1-5 - V*( - AX(2)/E(V)D - AX (2) 
z Py (z) Bay, {3) 


which is the PGF of the stationary queue size distribution due to server’s 
vacation period. It may be noted here that this result is obtained by 
Borthakur and Choudhury (1997). In particular, if we take limit of p — 0 
(i.e. no offered load in the system) in equation (3), we get 


lim Py(z) = {1 -V*(A~AX(z))]/E(V)IA - AX (2) 


= (Q(z) вау. (4) 
Hence the relationship between equation (3) and (4) is 
Py(z) = (1 — p)¢(z). 


Note that equation (4) represents the РСЕ of the number of customer that 


arrive during the residual vacation time for & multiple vacation of batch 


12 G. Choudhury 


arrival queue and it is independent of rate of service [for instance see Lee 
et. al. (1994)]. 
Similarly, taking limit of z — 1 in equation (2), we get 


lim P*(&.2) = (1 - р)[1 - V*(O]/6ECV)) 


which is the LST of the PDF of the length of a vacation period. 

It may be verified that the steady state probability that the server 
remains idle due to sever's vacation period is independent of the vacation 
time and equals (1 — р), which is the steady state probability of the server 
remains idle in a simple Мх /M/1 (or more general GI* /G/1) queue. With 


simple algebra, we have 
P,[The server remains idle] = lim P" (0; z) 


= (1-p) 


Let Ly, be the mean queue size during a vacation period when the system 
is idle, thén differentiating equation (3) with respect to z and then taking 
limit z — 1, we get on using L'Hospital's rule 


Ly 


dPy(z)/ а) = 
A(1 — p) E(X) E(Vn), 


where E(Vg) = E(V?)/2E(V) is the mean of residual vacation time. 
Also, it is well known to us [e.g. see Medhi (1997), gection-1.7] that 


. P.[System is idle] E(residual vacation time) 
(1 — p)E(V*)/2E(V) 
= W, (say), 


where W, = Expected time in the queue due to server’s vacation. 
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Thus we have 
№, = L,/AE(X) 


This shows that 


The mean waiting time due to server’s vacation 


(Mean queue size due to 
server’s vacation) 


(arrival rate) 


which justifies the Little’s law. 


Concluding remark 


Since vacation time is independent of the rate of service or more specif- 
ically service time distribution, therefore equation (4) holds good for even 
more general М” /G/1 types of queueing system. 
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Unsteady Rotational MHD Motion of a Dusty Viscous Fluid 
Contained in the Semi-infinite Circular Cylinder 
due to an Initially Applied Impulse on the Surface 


U.N. Das and A. Aziz 


Abstract: The unsteady rotational MHD motion of an electrically 
conducting dusty viscous fluid contained in the semi-infinite circular 
cylinder due to an initially applied impulse on the surface in the 
presence of a transverse magnetic field has been studied. Initially 
the liquid and dust particles are at rest. The motion is set up by an 
impulse of the shearing force on the surface. The governing equations 
Subject to the appropriate boundary and initial conditions have been 
solved analytically. The effect of mass concentration of dust particles 
on the fluid flow in presence of magnetic field and the effect of the. 
strength of magnetic field on velocity profiles at fixed time has been 
depicted graphically. 


1 Introduction 


Studies on the influence of dust particles on viscous fluid flows are of im- 
portance in petroleum industry and in the purification of crude oils. Other 
important applications involving dust particles in boundary layers include 
soil solvation by natural winds, lunar surface erosion by the exhaust of a 
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landing vehicle and dust entrainment in a cloud formed during a nuclear 
explosion. The unsteady motion of fluid resulting due to the pure rotation 
of a solid boundary or due to the application of & uniformly distributed 
shear stress along a solid boundary is of both theoretical and practical 


significance in fluid mechanics. 


À number of workers studied both steady and unsteady two dimen- 
sional axisymmetric rotational flow of a viscous fluid in view of its growing 
importance in viscous technical problems. Khamrui [1] analysed the slow 
steady motion of an infinite viscous flow due to rotation of a circular cylin- 
der. Iben [2] considered the non stationary, plane circular-symmetric flow 
of & viscous fluid which forms itself within as well as outside a rotating 
infinitely long circular cylinder. Bhattacharyya [3] studied the rotational 
motion produced in an enclosed field, contained in a circular cylinder of 
infinite depth. The disturbance was generated on the surface of the fluid 
by an impulsive couple. Later Mukherjee and Bhattacharyya [4] studied 
the rotational flow of viscous fluid due to the rotation of a circular cylinder 
or by the action of shearing stress on the boundary. Saffman [5] studied 
the stability of the laminar flow of a dusty gas with uniform distribution of 
dust particles. Michael [6] considered the Kelvin-Helmholtz instability of 
the dusty gas. Michael and Miller |7] has discussed the motion of dusty gas 
.enclosed in the semi-infinite space above a rigid plane boundary. Using the 
formulation of Saffman [5] many authors studied a number of gas problems 
and the results are well documented in a review by Marble [8]. Recently, 
Mandal, Mukherjee and Mukherjee [9] analysed the rotational motion of a 
dusty viscous fluid contained in the semi-infinite circular cylinder due to 
and initially applied impulse on the surface. In the present investigation 
we have extended the problem studied by Mandal et. al. [9] to the MHD 
case. The flow of the dusty gas is initiated by an impulsive shearing force 


n the surface in presence of a transverse magnetic field. 
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2 Formulation of the Problem 


Following Saffman [5] the basic equation to represent thé\motion of a dusty 
ffuid in presence of uniform magnetic field are given by 


p Ia" Gi voa] c -grad p* + ру? - KN(U —d)--Jx B 
(2.1) 
Ov В Aet 
m ж + бу) = K(v-u) (2.2) 
divi = 0 | (2.3) 
^ + ау Nð = 0 (2.4) 


where the velocities of the fluid and dust particles are and 7 respectively. 
N is the number density of dust particles each of mass m. Here K is 
the Stokes coefficient of resistance (for spherical particles of radius d, it ia 
67d), t, p*, p, ц, the time, the pressure, the density and the viscosity of the 
fluid respectively. The last term on the right hand side of (2.1) represent 
the force on the fluid due to the interaction of magnetic induction B and 
the electric current J in the fluid. 


In the present analysis, the following assumptions are made: 


1. The dust particles-are spherical in shape and are uniformly dis+ 
tributed. 


2. The interaction between particles themselves is not considered. 
3. The flow is fully developed. 


4. The number density of dust particle is constant throughout the mo- 


tion. 


5. The applied magnetic field is considered to be along the axial direc- 
tion of the cylinder. 


| 
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6. The flow induced magnetic field is neglected. 
7. There is no external applied electric field. 
8. Only the electromagnetic body forces are present. 


9. Fluid properties are invariable. 


Initially the fluid and dust particles are at rest. Consider the flow of 
a dusty fluid in a long circular cylinder of radius a. Disturbance is set up 
by an impulse of the shearing force on the surface. Referring the problem 
to cylindrical polar co-ordinate (т,0, х) we take the z-axis along the axis 
of the cylinder and the origin on the surface of the fluid. The prescribed 
impulsive shearing force will obviously generate a rotational velocity in 
the dusty fluid. Thus the fluid as well as dust particles will be moving in 
circular paths around the axis of the cylinder. Also a uniform magnetic 
field B(0, 0, Bo) is applied along the z-axis so that Bo will be transverse to 
the flow field, r denotes the radial direction measured outward from the 
z-axis. 


The symmetry consideration gives 


uj = из = 0, vy = v3 = 0, ug = ux(r, z,t) 


д 
v9 = vo(r, z, t) and 38 = 0 (2.5) 


where wy and vy are circumferential velocities of liquid and dust particles 

respectively. Since the distribution of dust particles is uniform, the number 

density N of the particles equals No, a constant throughout the motion. 
Then with the help of (2.5), equations (2.1) to (2.4) become 


дш _ „|Zu  1дш_ш ба 
a [92 тт mr 82 

+ KNo(v2 — u2) — c Віш (2.6) 

mt = K(ug 2 v2) (2.7) 
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where c is the electrical conductivity. 


Initial boundary conditions are 


мо = 0 at t < 0 for all z,v; = 0 att < 0 for all z (2.8) 


Ou 
Poz\2=0 = pL (22) » (2.9) 


is prescribed as function of r and t for т € a, t > 0, where Ро, is the 
shearing stress. 
цо — О as z — оо (2.10) 


из), а = 0, 2 20 (2.11) 


Introducing the following non-dimensional quantities и = (иза) /и, v = 
(v3a)/v, tq = t/t, тр = т (ит), аду = z/(vr)!?, ау = a/(vry?, f = 
(mNNo)/p (mass concentration of dust particles), т = m/K (relaxation ' 
time of dust particles) in (2.6) to (2.11), же get (dropping suffices) 

ди u 1дч u u 


_ GU, +0u u оч ua M2 

9p. ^ m ro v Wo у= 
(2.12) 
= = u-—wv (2.13) 
u=0 at t < О for all z,v =O at t < 0 for all z (2.14) 
Po:|z=0 = 5 = F'(r).6(t) (prescribed) (2.15) 

Oz Jz=0 

u- 0 as zoo (2.16) 
ul, = 0, 2 >0 (2.17) 


1/2 
where М = Bo (zz) is the Hartmann number and ó(t) is the Dirac's 
delta function. 


3 Method of Solution 


"We solve the present problem by using technique of Laplace transform. 
Taking Laplace transform of equations (2.12), (2.13) and using (2.14) we 
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get 
E au, 100 u Ou "a " 
pu = эз Кә лы damned (3.1) 


where ї = fy ue dt, = [у ve РФ, Re(p) > 0 
Eliminating 7 from (3.1) and (3.2) we get 

Ou lau u дт f 8 

OF la gta- [p( ti) мао (3.3) 


The solution of the above equation can be expressed in the form 


u(r, z p) = Ут. (3.4) 
n=l 
where u, = Л(олт)ф,(2,р) (3.5) 


and an can be determined form the equation (2.17) as the positive root of 
J (ana) = 0 (3.6) 


On substituting the value of &, from (3.5) in (3.3) we get 


To - [oe Meo (T +1)|%,=о (3.7) 


Solution of (3.7) subject to the Laplace transform of the condition (2.16) 
| 4 2 2 f He 
jm Аер -zf +M e»(1- ү) | (3.8) 
where A, is independent of z for all n. 
From (2.15) transformed-shearing stress on the surface is given by 


Poz\2=0 = F(r), T < a 


y CrJi(anr) (3.9) 
n=) 


p 23678 
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where 
Cy = 2a"? Шш? f (ж лд? (3.10) 
И 0 
Again from (3.4), (3.5) and (3.8) we get 
= OC f 1/2 
Pg, = — V Ji(anr) An lod + М2 +p (1+ г) (3.11) 
n=1 1 +р 


we now compare the expressions (3.9) and (3.11) to obtain 
2 2 f ps 
A, = —C, E + М + р (1 + i2) (3.12) 
The expression for velocity becomes 


= 1 2 2 f DS 
и = - 2, Cii (nr) +a *»(1e1L 


es[- [sore M ew 


where L^! is the inverse Laplace transformation operator. 


Inverse Laplace transform of (3.13) presents some difficulties and we 
restrict ourselves to calculate the velocity expression for large values of 
time t only. For large time t, p < 1 and 
2 2 Í um 2 2 1/2 
lod + M *»(1 1) ~ [o + М? + p(1 + f)" 
Then 


L E + M? +р (1 + o y 
: 1+р 


ехр E foz + М? +р (1 4 i5) | 


2 
= rum Бхр dmt | { (at) “We 0 nre 


22 UN Bascand A visis 


Equation (3.13) is then given by 


u = ea peer + DI 
x y СһЛ(оһт) exp | erem (3.14) 


This solution satisfies initial and boundary conditions given by (2.14) to 
(2.17) 


4 Particular Cases 


Case — 1: Motion due to impulsive shearing force applied within a circular 


area on the surface : 


We take 


Е(т) er for O0<r<b 


0 for b<r<a (4.1) 


Relation (4.1) corresponds to the situation where the applied force is 
acting within a circular area r = b, the rest of the surfaces being kept free 
from the impulse. 


From (3.10) we get 
C, = 2eb 77 (ала) Jo(anb)/(ana”) 


on substituting the value of C, in (3.14), we get the expression for velocity 
as 

—2eb? = 
aet + Pp > 


= -at 
x 2. Jo(cnb) J3 (osa) Ji (anton? | exp k a] (4.2) 


2(1 + f)/(4£)] - x |5 = z5] 


1+f 
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In absence of magnetic field (M = 0), the expression for velocity profile 
(equation (4.2)) becomes same as was deduced by Mandal [9] and in case 
when the fluid is clean (f — 0) and (M — 0) the expression (4.2) becomes 
same as was deduced by Bhattacharyya [3]. 


Case - 2: Flow due to applied impulsive force distributed over the 
circumference of the circle r = b; b « a: 


We take 
F(r) = 86(r — b) (4.3) 


where s is constant and 6 is Dirac delta function. 


From (3.10) we get 
Cn = 2a ?sbJi (Qnb) Jz (ana) (4.4) 


On substituting the value of C, in (3.14) we get the expression for 


velocity as 





bs _м? 
= MEET exp[-z^(1 + f)/(4t)] - exp E 
x ^» Л (олт) Л (оъ Б) Jy” (ana) exp (m т) (4.5) 


In absence of magnetic field (M = 0), the expression for velocity (equa- 
tion (4.5)) becomes same as that deduced by Mandal [9]. Again, when the 
fluid is clean and there is no magnetic field, that is, (f = 0) and (M =0), 
the expression (4.5) becomes same as was deduced by Bhattacharyya [3]. 
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0.0 02 0 4 0.6 08 10 
——— r 
FIG. 1: VELOCITY DISTRIBUTION OF DUSTY FLUID FOR DIFFERENT 
VALUES OF f WHEN t=), €»1,z-1,021, b 20 5,M 1,2, 


5 Discussion of the Results 


In order to illustrate the effects of mass concentration of dust particles on 
the flow field in presence of magnetic field and the effects of the transverse 
magnetic field on the velocity for two particular cases, numerical calcu- 
lations are carried out for suitable values of the parameters entering the 


problem and depicted in figures (1) and (2). From figures 1 and 2 it is ob- 
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served that magnitude of velocity of dusty fluid increases with the increase 
of mass concentration of dust particles for a fixed value of the Hartmann 
number M. From these figures it is also concluded that magnitude of ve- 
locity of dusty fluid decreases as the strength of the magnetic field increases 
for a fixed value of f. 





0.0 0.2 0.4 0.6 0.8 1.0 
— r 


FIG, 2: VELOCITY DISTRIBUTION OF DUSTY FLUID FOR DIFFERENT 
‘VALUES OF f WHEN t=1,5s=1,z=1,a=1, ba 0:5,M -1,M-2, 
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Weighted Composition Operators on С,(Х, E) 


Kamaljeet Kour and Bhopinder Singh 


On any space of functions with some structures there are two natural 
types of operators. one is the operator of multiplication and the other 
is the operator of composition. When we combine them somehow we get 
another type of operator known as the weighted composition operator. The 
operators belonging to this class has been the subject matter of systematic 
study for the last several decades on different function spaces, especially on 
ІР spaces, Н? spaces and spaces of continuous functions and have played 
a very significant role in the study of operators on Hilbert spaces. For a 
detailed information on such operators, we refer to the monograph [8] of 
singh and Manhas. 

This paper is & continuation of the earlier paper [7] in which one of 
the authors has studied composition operators on C,(.X, E). Here we shall 
report compact weighted composition operators on Cj( X, E). 


Preliminaries 


Throughout this paper we shall assume, unless stated otherwise, that X 
is a completely regular Hausdorff space and E a non-trivial Banach space 
(over K) with the norm ||. ||, where K is the field of reals (R) or complex 
numbers (C). Let С(Х, E) be the vector space of all continuous E valued 


?1991 AMS Subject Classification: 47B38, 46E40 
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functions on X, and let C,(X,E) be the subspace of C(X, E) consisting 
of those functions which are bounded. Then C4(X, E) is а Banach space 


under the norm 
[fll = sup{||f(x)|| : x € X). 

In case E = К, we simply write C,(X) instead of Cy(X, К). For any 
f in C(X) and any t in E, the function f, defined by setting f;(r) = 
f(z)t for all z in X clearly belongs to СЪ(Х, E). In particular, for f = 1, 
the constant one function on X, 1, € Съ(Х, E). By B,(E) (respectively, 
B,(E)) we denote the space B(E) of all operators (i.e. bounded linear 
transformations) from E into itself when it is equipped with the strong 


(respectively, uniform) operator topology. 
Weighted Composition Operators 


Let т be a B(E) valued function on X, Т a selfmap on X and let F(X, E) 
denote the space of all functions from X into E. Then the map f — v.foT 
is a linear transformation from С,(Х, E) to F(X, E), where the composite 
product «.f о Т is defined as 


r.f o T(z) = n(z)(f(T(z))) 


for all х € X. This linear transformation is denoted by И, т and is called 
the weighted composition transformation induced by the pair (7,1). In 
case W, т takes C (X, E) into itself and is also bounded, it is called a 
weighted composition operator on C;(X, E). | 

If T is the identity map on X, the corresponding operator W,7 on 
Съ(Х, E) is denoted by M, and is called the multiplication operator in- 
duced by т. On the other hand when т is the constant map and я(х) = I, 
the identity operator én E, for all т € X, W, т is denoted by Ст and is 
called the composition operator on СЪ(Х, E) induced by Т. For details, we 
refer to the recent survey article [9] of Singh and one of the authors. 
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Before characterizing those functions т: X — B(E) and T: X —^ X 
which induce weighted composition operators on СЪ(Х, E), we first note 
the following: 

If Wz p takes СЪ(Х, E) into itself and t € E, then W, т, is an element 
of Cy(X, E) and so we have 


| 


вир{|т(2)4|:=є X) вир{[[#у,гы (z)] : £ € X} 


|W ти lloc « oo 


Also, if x(X) is a bounded subset of B( E), we have 


1 = suptllm(z)]pus : x € X|| < oo. 
Further, if {£a} is a net in X such that ta — x in X, then we have 
Im(za)t — т() | = ПИ т (za) — W, тч (x)| — 0 


for all t € E. This implies that the map т: X — B(E) is continuous in 
the strong operator topology, that is, т € C( X, B,(E)). Note that 7 is not 
necessarily continuous in the uniform operator topology (see Jamison and 
Rajagopalan [4] for example). The continuity of т in the strong operator 
topology shows that the set N(x) = {т € X : x(x) Æ 0) is open in X. 
Also, we shall see that the map Т is continuous on N(x). But Т is not 
necessarily continuous on X X N (x) (see Jamison and Rajagopalan [4] again 


for example). 


Theorem 1 Let И, т be a weighted composition operator on C,(X, E). 
Then 


(i) x € C(X, B,(E)), 


(it) T is continuous on N (v). 
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Proof. Suppose that W, т is an operator on С(Х, E), where 7: X — 
B(E) and T : X — X. Then, as noted above, т € C(X, B,(E)). Further 
to prove the condition (ii), we assume on the contrary that there exists 
an ту in N(r) at which Т is not continuous. Then there exists a net 
{xq : а € A} in N(x) such that х, — zo but T(xq) does not converge 
to T(z). This further implies that there exists an open neighbourhood G 
of T'(zo) such that, for each ag € A, T(za) @ С for some a > œp. Thus 
we can find a subnet {yg} of {xa} such that T(yg) Є С. The complete 
regularity of X gives an f € C,(X) such that 0 € f <1, f(T(zo)) = 1 and 
f(X\G) = (0). Let t be a non-zero vector in E such that v(zo)t 7 0. and 
consider the function f; € СЪ(Х, E) (as defined earlier). But for each D, 


we have 


|| т. (ув) — Ws felso) 
= ||л(ув) (7(у8)) — (20) £f (T(zo))ll 
: = [mn(zo)t|| > 0. 


This implies that W,.rf; ¢ Ca(X, Е) which contradicts the fact that Wr 


is an operator on C,(X,E). Therefore, Т must be continuous on N (v). 
Remark. Conditions (i) and (it) in Theorem 1 above are necessary but 
not sufficient for Их т to be a weighted composition operator on C,(X, E) 


as the following example shows: 


Example. Let X = R, E = L?(R) (with Lebesgue measure), 
l if 0 


and a(x) = 7, — I, where т, — Г is the operator of translation by z. Define 
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f in СХ, E) as follows: 


so | V Iz[xc. si) if |v| » 1. 


X(z,74-1) if 0< |x| <1 


Then note that И, rf(0) = 0 and IW. rf(i)la- V2, n = 1,2,.... This 


implies that W, т is not an into map. 


Remark. For sufficiency, we must assume that т(Х) is bounded and T 


be continuous on N (v). 


Compact Weighted Composition Operators 


Recently many authors have studied the compact weighted composi- 
tion operators on continuous function spaces. Kamowitz [5] characterizes 
the compactness of these operators on the Banach algebra C(X) of all 
complex valued continuos functions on a compact Hausdorff space X. His 
result has been extended to a larger class of spaces of continuous functions 
by Singh and Summers [10], and Jamison and Rajagopalan [4]. Further, 
Chan [1] and Takagi [11] have respectively improved upon the result of 
paper [4]. Feldman [2] deals with compact weighted composition opera- 
tors on Banach lattices while Lindstróm and Llavona [6], and Gonzalez 
and Gutierrez [3] describe the relationship between compact and weakly 
compact homomorphisms (of composition type) between some algebras of 
continuous-functions and differentiable functions, respectively. . 

Following Theorem 2 of J amison and Rajagopalan [4], we shall present 
in the following theorem a characterization for compact weighted compo- 
sition operators on Cj( X, E). This theorem, in turn, extends Theorem 2 
of Feldman [2] and à Theorem of Singh and Summers [10]. 
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— Theorem 2 Let Wk r be a weighted composition operator on СЪ(Х, E). 


Then И т ts compact implies that the following conditions hold: 


(2.1) п is continuous ın the uniform operator topology; 
(2.2) Each x(x) is a compact operator оп E; and 


(2.3) T(N(m,c)) is finite for every є > 0 
where №ъ№т,є) = (x € X : ||n(z)|| > є). 


Proof. Suppose that W, is compact. Then the proof of condition 
(2.1) is the same as in the proof of Theorem 2 of [4]. To see that (2.2) 
holds, let хо € X be fixed and let {tn} be an arbitrary sequence in 
Eo —' [t € Е: |tt| € 1). For each n, consider the constant t,-function 
1,, (as defined earlier) in Сь(Х, E). Then (1, : n > 1) is a bounded se- 
quence in C)(X, E), and the compactness of Wy т implies that {Wren } 
has a convergent subsequence in Cj(X, E). But И Tin (20) = T(zo)ta. 
Therefore we conclude that 7(zg)t, also has a convergent subsequence in 
E. Thus v(xz9) is a compact operator on E. 

Next, assume that (2.3) is not satisfied. Then there exists an є > 0 
and a sequence {a;} of distinct points in T(N(z,e)). As constructed in 
the proof of Theorem 1 of Feldman [2] we choose a collection of distinct 
neighbourhoods {N,} for points {a;,} in {a;}, and choose a collection of 
continuous functions {f,} such that 0 < f, € 1, fi(a,) = 1 and f,(X \ 
№) = {0}. For non-zero vectors t, to be chosen presently in E, define the 
function h; as h;(z) = f,(x)t, for alla € X andj > 1. Then (^, : j > 1} is 
a bounded sequence іп C (X, E), and so by compactness of Ў.т, {W т^) 
should have a convergent subsequence. But for т € X and ё, so chosen 


that [|m(z)t,|| > є and T(z) = a,, we have 


IV 


Ir), (T (z)) — (а) (T (a) 


= |Ir(z)tll-1(a,)1 > € 


ПИ r = Wr M 
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a contradiction. 


Remark. 'The condition of Theorem 2 above are necessary but not suff- 


cient as can be seen by the example given below: 


Example. Let X = {0,1,2,...} with the discrete topology and let E 
be a separable infinite dimensional Hilbert space with ep, e1,é€2,... as an 
orthonormal basis. Let T': X — X be given by T(z) = 0 for all z c X; 
and for each x in X, define л(т)» = (v,€r)ez for all v € E. Then the 
varification of conditions (2.1)-(2.3) are trivial. But И» т is not compact 
because if we define f, : X > E as fi(r) = ej for all x € X, where 
t=(,1,2,...; then (fi) із a bounded sequence in C,(X, E) with || fill = 1 
and И, тў —W2,T fill =1 fori 4 j. 


References 


[1] J. T. Chan, Operators with the disjoint support property, J. Oper- 
ator Theory 24 (1990), 383-391. 


[2] W. Feldman, Compact weighted composition operators on Banach 
, lattice, Proc. Amer. Math. Soc. 108 (1990) 95-99. 


[3] M. Gonzalez and J. M. Gutiérrez, Compact and weakly compact 
homomorphisms between algebras of differentiable functions, Proc. 
Amer. Math. Soc. 115 (1992), 1025-1029. 


[4] J. E. Jamison and Rajagopalan, Weighted composition operator on 
C(X, E), J. Operator Theory 19 (1998), 307-317. 


[5] H. Kamowitz, Compact weighted endomorphisms of C(X), Proc 
Amer. Math. Soc. 83 (1981), 517-521. 


34 K. Kour and B. Singh 


[6] М. Lindstrom and J. Llavona, Compact and weakly compact homor- 
phisms between algebras of continuous functions, J. Math. Anal. 
Appl. 166 (1992), 325-330. 


[7] Bhopinder Singh, Composition operators on Cy(X, E), Indian J. 
Pure Appl. Math. 26 (1995), 97-101. 


[8] В. К. Singh and J. S. Manhas, Composition operators on function 
spaces. North Holland Math Studies, 179, Amsterdam, 1993. 


[9] В. К. Singh and Bhopinder Singh, Compact weighted composition 
operators on spaces of continuous functions: A Survey, Ezxtracta 
Math. 10 (1995), 1- 20. 


[10] R. K. Singh and W. H. Summers, Compact and weakly compact com- 
position operators on spaces of vector valued continuous functions, 
Proc. Amer. Math. Soc. 99 (1987), 667—670. 


[11] H. Takagi, Compact weighted composition operators on certain sub- 
spaces of C(X, E), Tokyo J. Math. 14 (1991), 121-127. 


Dept. of Mathematics Dept. of Mathematics 
University of Jammu Govt. College of €ngg. and Technology 
Jammu ~ IBO 004 Jammu - I8O ОО! 


The Bulletin, GUMA 
Vol. 4 (1997) 85-41 
Published: August 1998. 


Generalized Monotone Operators and 
Complementary Problem 


Sudarshan Nanda 


Abstract: The concept of F-monotonicity was first introduced by 
Kato and this generalizes the notion of monotonicity introduced by 
Minty. The purpose of the present paper is to define various types of 
F-monotonicities, discuss relationship among them and study their 


applications to variational inequalities and complementary problem. 


Keywords: Variational inequality, complementary problem, closed convex cona: 


monotone operators, coercive operators. 


1 Introduction 


Let X and Y be any two (real or complex) Banach spaces. Let F and T 
be two non-linear operators such that 


F : D(F) = X > R(F) CY 


and 
T:D(T) c X 5 R(T) сҮ", 


°This paper was written while the author was visiting International Centre for Theo- 
retical Physics, Trieste, Italy under the Associateship scheme. 
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whereD and R denote the domain and range of the operators, respectively. 
Ihe concept of F-monotonicity of the operator T' was introduced and dis- 
cussed in Kato [2,3]. The usual concept of monotonicity introduced by 
Minty [5] and discussed by many others is a special case of this concept, 
the case for which F is the identity operator. 

The purpose of this paper is to introduce various types of F-monoto- 
nicities and discuss the relationships among them. We first introduce the 
following definitions: 


T' is said to be F-monotone if 


> 


Ке(Тт — Ty, F(z — y)) > 0 for all x,y € D(T) 
and strictly PF-monotone if 
Re(T'z — Ty, F(x — у)) > 0 for all x,y € D(T),z zy 
T is strongly йоне if there exists а constant c » 0 such that 
Re(Tz — Ty, F(x — y)) > e||F(z — у)? 


for all x,y € D(T). Т is called a-F monotone if there exists a continuous 
strictly increasing function a : [0, оо) — [0, оо] with a(0) = 0 and a(r) — 


oo as r — oo such that 
Re(T'z — Ty, F(x — y)) > c||F(z — y)le(lF(x — y)||) 
for all z, y € D(T). 


T' 18 F-coercive if 


(Tx, Ty) 
|| Fl 


Let Y = X and F be the identity operator on X. Then the above 


concepts are just equivalent to monotonicity, strict monotonicity, strong 


— oo as ||F'x|| — oo. 


monotonicity, a-monotonicity and coercivity in the usual sense. 
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Let X bereflexive, Y — X* and X* be strictly convex. Let F : X — x* 
be the duality map in the usual sense (which is well-defined if we assume 
Fæ] = |||). In this case F-monotonicity of T (which has range space 
Y" = X) means that Т is accretive in the sense of Browder [1]. It may be 


noted here that in this case F is itself monotone in the usual sense. 
Let K be a convex subset of X. Then 


r€K:(Tz,F(y-z)»0foralye K (1) 


will be called a generalized variational inequality and any z € K which 
satisfies (1) will be called a solution of (1). Similarly 


z EK : (Ту, F(z — у)) > 0 for all y € K (2) 


will be called another generalized variational inequality. Note that if F is 
the identity map in (1) and (2), then we get the usual variational inequal- 
ities. Let 5; and So, respectively, denote the sets of solutions of (1) and 
(2). | 
Let K bea convex cone in X. Let K* be a subset of Y* defined by 


K^ = {y € Y* :(y, Fx) > 0 for all y € К}. 


Then 
rz € К, Tz e K*, (Tz, Ех) =0 (3) 


will be called a generalized complementary problem and any т € К which 
satisfies (3) will be called a solution of (3). Note that if F is the identity 
map, then (3) reduces to the usual complementary problem. Let С denote 
the set of all solutions of (3). 

Before describing the results we first quote the following definitions 
which will be needed in the sequel. The mapping F is said to be positive 
homogeneous if 

F(tz) —tF(zx)fort»0 
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Е is called symmetric if F(x) = F(—z) and antisymmetric if F(x) = 
—F(-z) for all x € D(F). F is said to be additive if 


Е(ж + у) F(z) + F(y) forz,y € D(F). 


The mapping Т is said to be hemicontinuous if (Т) is convex and for any 
x,y € D(T), the map t > T(ta + (1 — t)y) of [0,1] to Y* is continuous for 
the natural topology of [0, 1] and the. weak topology of Y*. 


2 Results 


Theorem 1 ГТ: D(T) C X — Y* is a-F monotone, then it is 
strictly F-monotone (hence F-monotone) and F-coercive. In par- 
ticular, every strongly F-monotone operator ts strictly F-monotone 


and F-coercive. 


Proof. It follows from the definition that every a-F' monotone operator is 
strictly F-monotone. To prove F-coercivity let Т be a-F-monotone. Then 
there exists a continuous strictly increasing function a : [0,оо) — [0, oo] 


with a(0) = 0 and a(r) — oo as r — oo such that 
Re(Tz — ТО, Fr) > ||Fz||o(|Faz|]) 
Or 
Re(T'z, Ет) > Re(T0, Fr) + ||Fzl|o(||Fz||) 
Or 


Re(Tx, Ет) _ Re(T0, Fx) 


Ste ee E 
Ра Е | 


— oo as ||Fz|| — oo. 


Therefore T' is F-coercive and this completes the proof. 
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Theorem 2 If Е is antisymmetric and Т is strictly F-monotone, then 


Sı is either empty or a singleton. 
Proof. Assume, to the contrary, that 2, and ze are in $1. Then 
re(Tz1,F(z;9 5271) 250 (4) 


and 
re(T'z5, F(x, — хз) > 0 (5). 


since F' is antisymmetric it follows from (4) that 
re(Tz;, F(x — 29) < 0 (6) 

From (5) and (6) we obtain 
re(T'z, — Tx, F(x, — z3) < 0 (7) 


Since T' is strictly F-monotone, this is impossible unless z, = т), and this 


completes the proof of the theorem. 


Theorem 3 Let F be positive homogeneous, let T be F-monotone and 
hemicontinuous. Then 
91 = So. 


Proof. Let x € Sj. Since Т is F-monotone 
re(Ty, F(y — x) > re(Tz, F(y — х) > 0. 


Hence т € So. 
Let x € S2. Let y € К be arbitrary. Since К is convex, 


и = (1—t)et+ty, 0ct«1 


is in K. Hence 


40 S. Nanda 


(Ту, F(yy — x) > 0. 


But y: — x = t(y — x) and since F is positive homogeneous we have 
(Ty, Fly m z) 20 
Since T' is hemicontinuous we now have 


(Tz,F(y-2z)20 


— 


and hence z € Sı. Thus Sı = 55 and this completes the proof. 


Theorem 4 Let K be a closed conver cone in X, F additive and antt- 


symmetric. Then C = 5). 


Proof. Let т € Sı. Таке y = 0. Then (Tz,F(—z)) > 0. Since F is. 
antisymmetric we have (T'z, F(x) < 0. Take y = 2z € K. Then (Tz, Fx) > 
0. Thus (Tz, Fx) = 0. Assume, to the contrary, that Tx є K*. Then there 
exists yo Є К such that (Tx, Куо) < 0. Since т € Sı and F is additive we 
have 

0 > (Tz, Fyo) > (Tz, Fx) = 0. 


This contradiction proves that Tz € K* and hence т € C. Thus 5] C C. 
Since F is additive, it is obvious that C C 5. Thus C = S; and this 


completes the proof. 
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Hydromagnetic Divergent Channel Flow of a 
Newtonian Electrically Conducting Fluid 


Deva Kanta Phukan 


Abstract: The two-dimensional boundary layer flow through a 

divergent channel of a Newtonian electrically conducting fluid in the 

presence of transverse magnetic field is studied. A special form of 

magnetic field is chosen so as to yield similarity equation. First, 

perturbation solutions for the case of large magnetic parameter are 

derived. Finally, a direct numerical solution of the resulting boundary 

value problem is obtained using a Runge-Kutta shooting algorithm 
with Newton iteration. | 


Keywords: Divergent Channel, Electrically Conducting fluid, Mag- 
netic Field. 


1 Introduction 


The divergent flow problem between two non-parallel planes was solved by 
Jeffery [1] by reducing the problem to an elliptic integral equation. Sri- 
vastava [2] extended this problem to an electrically conducting fluid in the 


presence of transverse magnetic field. He found that with the application 
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of magnetic field it is possible to have purely divergent flow without any 
secondary flow for greater angle between the two planes. Rosenhead [3] 
presented the solution of two-dimensional incompressible laminar flow in a 
diverging channel with impermeable wall. Terril [4] forwarded the solution 
obtained for small laminar flow through the divergent channel with suc- 
tion at one wall and blowing at the other wall. Srivastava and Hazarika [5] 
discussed the problem [2] bý, shooting method, | 

A solution for the flow through a convergent channel which is based 
on the boundary layer (approximation) equation was first obtained by 
Pohlhausen [6]. The solution was performed in a cartesian Co-ordinate sys- 
tem hence there are two non-zero velocity components along two mutually 
parpendicular directions. Goldstein [7] presented an interesting discussion 
of this problem. 

The present work deals with the flow of an electrically conducting fluid 
through a divergent channel which is based on boundary layer approxi- 
mation in the presence of transverse magnetic field. Similarity solutions 
are possible only if the magnetic field takes a special form. Linearized so- 
lutions for the case of large magnetic parameters are presented. Finally, 
a direct numerical solution of the similarity equation is obtained using a 
Runge-Kutta Shooting algorithm with Newton iteration. 


2 -Basic Equation 


The basic boundary layer equations for the steady two-dimensional flow of 
an electrically conducting, viscous, incompressible fluid in the presence of 
а magnetic field B(x) are taken as (see, Djukic [8]); for Newtonian fluid) 


Ou ду 


ана 2.1 
Oz ' By 0, (2.1) 
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du OU да  cB*(a) 
—— — —Q NER S 2.2 
| ug, t UA z BAT (U — ч), (2.2) 
subject to the boundary conditions 
у = 0: u = 0,0 = 0; у = оо: u = U (g), (2.3) 


where u and v are the flow velocities in the z— and y— directions respec- 
tively: U is the velocity component outside the boundary layer, p the fluid 
density. v the kinematic viscosity and c the electrical conductivity. In 
equation (2.2), the secondary effects of magnetic induction are ignored, i.e. 
the induced magnetic field is neglected being small in comparison to the 
applied magnetic field. Furthermore, we assume that the external electric 
field is zero and the electric field due to polarization of charges is also 
negligible. 
As in Sinha and Choudhury [9]. the potential flow near the source is 
taken to be 
(а) = 2, (2.4) 
T 


which with uj > 0 represent two dimensional divergent channel flow and 
leads to similar solution. 
We introduce the following change of variables (see, Schlichting [10]) 


n(z,y) = ym = ae (2.5) 
V(z,y) = J/U(z)vzF(n) = /vuiF() (2.6) 


We obtain the velocity components 


uc а = U(z)F'(), v = -2E = yrn Ра) (2.7) 


It can be easily verified that the continuity equation (2.1) is identically 
satisfied. Similarity solution exis:s if the magnetic ficld has the special 
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form (see, Chiam [11]) 
B(z) — A (2.8) 


Using (2.4)—(2.8) in (2.2) gives the following equation 

F" 4 F^ - М(1- F')—1-0 (2.9) 
subject to the boundary conditions | 

F(0) = 0, F'(0) = 0, (оо) = 1 (2.10) 


where М = o BZ /puy, is the hydromagnetic parameter. 


3 Approximate Solution for Large M 


We may obtain a linearized solution by considering the magnetic parameter 
M to be large. The procedure given below essentially follows Sapunkov's 
work [12] which was concerned with studying the problem of an electrically 
conducting non-Newtonian power-law fluid in the presence of transverse 
magnetic field which has recently been followed by Chiam [11] to present 
the Perturbation solution of hydromagnetic flow over a surface stretching 
with a power-law velocity distribution. We first assume 


z = V Mm, f(z) = V MF() (3.1) 
which implies ў dF 


Equation (2.9) then becomes 
f"(z) — f(z) +1 = e — f^(z)) (3.3) 
where є = 1/M. The corresponding boundary conditions are 


/(0)-= 0, f'(0) = 0, f'(oo) = 1. (3.4) 
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The unknown function f(z) is expanded in terms of powers of the small 


parameter є as follows: 


f(z) = folz) + efi(z) +? fa(z) +... (3.5) 

Equation (3.3) with the help of (3.5), then gives the following set of equar 
tions: 

fo -fo = -1, (3.6) 

Л = 1- fy’, (3.7) 

2—-/3 = 2р. (3.8) 


The above uncoupled linear equations can be readily solved in sequence. 
The appropriate boundary conditions are 


fo(0) = 0, fo(0) = 0, /0(со) = 1, (3.9) 
f1(0) = 0, fi(0) = 0, Л(оо) = 0, (3.10) 
f2(0) = 0, f2(0) = 0, f2(oo) = 0. (3.11) 


The solution of equations (3.6)-(3.8) satisfying the respective boundary 
conditions (3.9)-(3.11) are 


fo(z) = z-Fexp(-z)- 1, (3.12) 
h(z) = (04 + 6z) exp(—2) + exp(—2z) — 5}, (3.13) 
falz) = ag { (182? + 42z + 23) exp( —z) 

+ 12z exp(—2z) + exp(—3z) — 38}. (3.14) 


From these solutions, we can easily compute the approximate skin friction 


co-efficient which, in this case, is given by 
2 2 
"(Оу=1——є— Ĉe 
f (0) a ae е, (3.15) 


to second order in the small parameter є. 
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4 Numerical Solution 


A numerical solution of equation (2.9) with the boundary conditions (2.10) 
cannot be obtained because to integrate a third order ordinary differential 
equation by the standard numerical methods, three conditions must be 
prescribed at the initial point of integration. Since one condition is always 
missing at this point, ie. the value of F"(0) is not prescribed at the 
initial point of integration, so it is necessary to determine the value of 
the missing initial conditions by an appropriate numerical method. The 
shooting method [13.14] is & suitable method which enables to convert 
a boundary value problem to an initial value problem. Hence a direct 
numerical solution of equation (2.9) subject to boundary conditions (2.10) 
can be obtained by a shooting method using the Runge-Kutta algorithm 
with Newton iteration in double precision. The trouble of this procedure 
is that we have to make an initial guess for the value of F"(0) to initiate 
the shooting«process. The success of the procedure depends very much 
on how good the guess is. The process i8 frequently very sensitive to 
the starting value. This problem is made worse by the presence of the 
ТРЕ А parameter М in the governing BOE In order to study the 
cases corresponding to many values of the parameter M, we would have to 
. make a great number of guesses. Fortunately, we have computed the value 


‘of Е" for various values of the magnetic parameter M. 


5 Results and Discussion 


The variations of dimensionless velocity F” with 7 across the boundary layer 
for different values of the magnetic parameter M, are shown in Fig. 1. It 
can be seen from Fig. 1 that Р” decreases with the increase of the magnetic 


parameter M. As M increases, the Lorentz force which opposes the flow 
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Fig. 1. The dimensionless velocity distribution F' for 
different values of magnetic parameter M. 


also increases and leads to enhance deceleration of the flow. That is, the 
effects of the magnetic field are to decrease the velocity of the а 
fluid across the boundary layer. 
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Fundamentals Concerning Water Waves: 
Part III : Solitary Waves 


Matiur Rahman 


[ Introduction: 

The contents of this paper is in continuation of the earlier two papers by 
the author published in Volume 2 (December 1995) and Volume 3 (December 
1996) with the titles Fundamentals Concerning Water Waves: Part I: Airy 
Waves and Fundamentals Concerning Water Waves: Part II: Stokes Waves 
respectively. In the first paper the author investigated the surface gravity waves. 
There he considered the liquid (water) as inviscid and the motion irrotational 
so that a velocity potential exists. Under these assumptions the laws of motion 
in water for three dimensional motion were derived by taking into account both 
the horizontal and vertical accelerations. Several interesting results concerning 
surface gravity waves, progressive waves, surface waves and group velocity were 
established. 

In the second paper. the author investigated various aspects of stokes waves. 
'This is essentially a premise of the non-linear wave theory or finite amplitude waves 
as a departure from the small amplitude (linear theory) wave motion discussed in 
the first article. Here also the motion is considered as irrotational but the velocity 
potential and the surface elevation 7 are expanded as series in terms of a suitable 
perturbation parameter є (= ka which is related to the ratio of the wave height 
to the wavelength). Results have been obtained upto second order perturbation 
quantities. Finally, the stability aspects of the water waves have been shortly 
discussed. The present paper has been devoted to solitary waves pertaining to - 
the shallow water waves. This requires dealing with non-linear theory of water 
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waves. The author has presented a mathematical analysis for two dimensional 
case of shallow water and constant depth wave propagation. 


– U.N. Das] 


1 Nonlinear long waves in shallow water 


In the Stokes second-order theory, the surface profile for small parameter 
kA (the wave steepness) has been obtained in (1:35). In very shallow water 
this Stokes wave profile becomes (see also Struik (1926)). 


3A*k 
n(x. t) = Asin(kz — ot) — 41) 952 8 egt) — —. (1) 


when kh — 0. We have scen in the previous section that for deep water 


waves when kh — oo this surface profile is 


| 2 
n(r,t) = Asin(kz — ot) — — cos 2(kz — ot). (2) 
- From this analysis it is clear that the nonlinear theory of the Stokes 
expansion for КА < 1 is valid for deep water waves only. However, from 
this nonlinear theory in the case of shallow water waves, we have i «1 
such that (kA) < (4) (kh)*. This subsequently yields 
A 


4 2 


Thus for this case both kh << 1 and 1 «« 1 must hold. Let us define 
these two parameters as 


y-—-kh««1 and TI (3) 


Here the second restriction is a much more severe one for many prac- 
tical problems. Therefore, it is essential to investigate a nonlinear theory 


of shallow water waves. Historically, two different theories were developed, 
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one by Airy (1845) and the second by Boussinesq (1877), which led to 
opposite conclusions. However, these differences were resolved by a fun- 
damental paper by Ursell (1953). In particular, Ursell (1953) showed that 


the ratio, 


_ Afh h kA 
M 4 
Un = 7 y 7 (y " 
which is known as the Ursell number, plays a very important role in decid- 


ing the choice of approximation. 


2 Mathematical developments 


A formal procedure for systematically developing shallow water wave the- 
ories was first introduced by Friedrichs (1948) and has also been outlined 
by Stoker (1957). The following mathematical equations for the constant 
depth case are mainly due to Benney (1966) and Peregrine (1967) and the 
analysis is made for the two-dimensional case. Because of the presence 
of two parameters, it will be convenient to express the equations in non- 
dimensional form prior to the perturbation procedure. For convenience we 
define the governing equations with their boundary conditions with up- 
per case symbols for the dimensional case and lower case symbols for the 
non-dimensional case. 


Laplace’s equation: 


9?Ф ‚9% 
ax? t azz 7° (5) 
_Dynamic surface boundary condition: 
04 0 06? _ 
ar ta ix) + (33) =° li (6) 
Kinematic surface boundary condition: 
0j 09605 Ae 


ƏT axax oz * 47-8 Л 
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Bottom boundary condition 
-5 =0 at Z=-h (8) 


Here Ф = velocity potential, ñ= wave elevation, X, Z = Cartesian coordi- 
nates, Т' = time, and g = acceleration due to gravity. | 

The linear theory suggests that the following dimensionless ‘variables 
are appropriate to make the above equations dimensionless. 


т=ЕХ, z= Z, t = k\/ghT 
| (9) 
i 6 ahs 
The velocity components are normalized as 
U = $ = A Van = {һи 
(10) 
W = $8 = (д) ав = (4) (4) va 
It is to be noted that the difference in scaling for horizontal and verti- 
cal components is required by continuity. The normalized equation and 


boundary conditions are: 


Laplace’s equation: 


p | 0$ 
aa 5: = 0, —1 < z < 7 (11) 


Dynamic boundary condition: 


(Seen) +5 (E + (32) | -o z=6n (12) 


Kinematic boundary condition: 


Y Fax +5{ ae =<, z= ôn (13) 
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-Bottom boundary condition: 
— —0 at z=-1 (14) 


~ 


Keeping 6 "T" we first-ássume that y = kh << 1. Since ф is analytic 
we-can expand it as a power series in the vertical coordinate so that 


(в, zit) = S (z + D" de (,t) (18) 
n=0 


Then we have 


д = д 
дф > (z + "22 


дт n=0 
EI = е + E 
20 a Ут + 1)%—10, 
= Уы + 1)" 18, 
= Ne + 1)(z + "6444 
E = Y n(n + 1) (2 + 1/7 ss; 
= E n(n + 1)(z + 1)"714, 


$n 1)(n  2)(z + 15,45 
n=0 с 


Using these equations in Laplace's equation (1t), we obtain 





PY + THs Vint а ача =0 (10) 
n=0 7 n=0 


“a 
a 
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Since 2 is arbitrary, —1 < z < бт the coefficient of each power of (z + 1) 


must vanish yielding a recurrence relation 


P bp 


Pea "(n * 102) Oa? 


n = 0,1,2,3, ... (17) 


But the bottom boundary condition implies that se = 0 at z = –1. 

That is У (п + 1)(2 + 1)%ф, 1 = 0 at z = —1. Corresponding to 
n = 0,4; = 0. It is clear from (17) that $1 = фз = $2441 = 0. Also from 
(17) , we have 


= 7 9% (18) 
п=4: фо 55H = Te 


If we assume that фу = O(1) then it may be calculated that фә = O(4^), 
фа = O(7*), ds = O(49), and so on. Thus with an error of O(y°), the 
velocity potential is 


b= фо Toni. (+ „озу a9 


Now we will investigate the dynamic and kinematic free-surface boundary 
conditions. Let us put б = 1 + 67. With this substitution (12) reduces to 


r [se Fas (ic) v 
(5) ve (8 59) 


51 
+ ce = 009) (20) 
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. Equation (13) can be written as 


dE + oe Ode -те е} 


60t дг \ ðr 2° 923 

po түй 391% 
— Aero Topsy YO 6 , 
= Yea tas әм tO) (21) 


Now define uy — Sin and then differentiating (20) with respect to x we 
obtain the dynamic boundary condition as 


b йш e T4 ty ы (52) 


Ot Or 802 
б 
-Eeu y ке (58) |= От) (22) 


The kinematic boundary condition (13) can be simplified as 


19$ | a Y^ „a8 uo 
6 at + Өт (s TT. s) | 
дщ P з  (Ou\ 4 

| xe - EA (54) = 005 (23) 

Once 7 and up are solved, the actual velocity components are given by 

Op — ү? 2 Fup 4 

- pros Au OQ) (24) 
ш = Ul —Ү (2+ е ~ + OQ) (25) 


The pressure field is obtained from Bernoulli's equation whose dimension- 
less form is exactly 


0d 909? 1 (8¢\? 
ZEUFE "e +a e) (26) 
. where p has been normalized as p = -^r. The approximate pressure field 
may be obtained by using (24) and б into (26) 


m дф y! 2 O ug 
—p = ZEE. 8t — 5 (2 +1) dxdt 
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-= ; 


+2 sfa- Y (z t 1)*uo 
+ 72(2+ 1)? (2) Y + O(^) (27) 
Equation (20) can be used to eliminate 9 so that | 
(бп — z) — Tie — (z+1)"} p 
+6 m = ($2) они (28) 


Now we can introduce the depth averaged horizontal velocity defined 
by 


<u> = l ” (#2) ae 


e? 
= _ 2 4 
= ug ¢ o2 + O(7*) (29) 
which can be easily inverted to yield 


ug =< u > Paa Lad < u » +О(4%) (30) 
6^ дт? 
Substituting equation (30) into eqn (23) it follows that (kinematic bound- 


ary condition) 
д6 
8t 


which is just the depth averaged continuity equation. 


+62. <u>) = 0(ү). (31) 


Now expressing up in terms of < u >, eqn (22) (dynamic boundary 
condition) and (28) (the pressure) can be written respectively, as 
Ə<u> д 18€ xe 


теш; кошы. рр ссср 
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im (s <u>) ] = O(4^) (32) 


<u> 
dxdt 


46 <u> P ке [эсш 
Өх? дт 


° 2 
р = (ón-2)- Zl? - (z+ 1)3] 





-O(**) (33) 





These equations are valid for 6 only. Extension to 6? is straightforward 
but tedious. Expressing equations (31), (32) and (33) in physical variables, 
we obtain 


a) ze) « U >) =0 (34) 


Ó д 9 - 

8T <U>+<U> a <U> t95y(h ii 
10 з. OF 

age (b+ a pgs <>) 


д ( 1 : д? 
*t3x (5 (h +ñ) (<U >) (зя < 0 >) 


ne 7)? (ax «U >) 
2 К \ ex 
1 


; (i 0-5, (эт <U >)) = 0 (35) 


апа 


-— 


P = pol 2) - [nes - ewe <u> 


9? 8 = ү? 
+ (<u> эж <0>2-(0.0))] (36) 


where 


A 0 д 
<U >= yyh < ч >, у = Ко | (37) 
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Limiting cases: 
(A) Airy’s theory for very long waves: y — 0 and ё = O(1) 
Airy's theory gives the leading-order approximation for very long waves of 
finite amplitude, because for very long waves y = kh — 0 and 6 = i ES 


O(1). Thus omitting terms proportional to гу? for equations (31), (32) and 
(33), we obtain in physical variables 


д И д - 

ap t 9) + 3x (0 +?) <U>}=0 (38) 
8 8 д Е 

ap <U> + <U> у <U> tag (ht) = 0 (39) 
P = рд(й— 2) (40) 


These equations are valid for variable depth h(x) also. A distinguishing 
feature of Airy’s approximation is that the pressure is hydrostatic. 
(B)Boussinesq theory: 6 = O(4?) < 1 


| For weakly nonlinear and moderately long waves in shallow water 6 = A < 
1 and +? = kh < 1, eqns (31), (32) and (33) reduce to 








д 

Ez 60) <u>} =0 (41) 
д<ч> д On 

8t docu UU pa | 

py 9? 

-balg <>) =0 (42) 

2 
Е Е y? o0 T 
= (ên — 2) + 162? +22) = а (48) 
In physical variables these can be written as 

LL NS 2G +h) <U>}=0 (44) 
oT’ ax” 
8<U> 8 99 


ату 


8T 8X ӘХ | 
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h? & (0 
a S P2 ð (OcU» 
Р = pg(ñ - 2) + £(2? + 2hZ) => T } (46) 


Equations (41), (42) and (43) or equivalently (44), (45) апа (46) are called 
the Boussinesq equations. Here the pressure term із по longer hydrostatic. 
Moreover, the Áiry and Boussinesq theories differ by the linear term mul- 
tiplied by +” in eqn (43). It is to be noted here that the two parameters 
6 = (2) and у = kh can be treated as nonlinearity and dispersive parame- 
ters respectively. The Boussinesq theory accounts for the effects of nonlin- 
earity 6 and dispersion ^^ to the leading order. However, when 6 >> 4?, 
they reduce to Airy’s equations which are valid for all 6; when 6 << 7? they 
reduce to the linearized approximation with weak dispersion. When 6 — 0 
and у? — 0, the classical linearized theory can be recovered. The foregoing 
analysis is valid for the constant depth case. For the variable depth case, 
this analysis can be extended easily (see Mei (1989)). Extensive work on 
the shallow water waves and the usefulness of the characteristic method 
has been reported by Rahman (1994) in his excellent book on water waves 
published by Oxford University Press and thus will not be repeated here. 


3 Nonlinear dispersive waves 


Mathematical expressions of the form exp{i(ka — ot)} = exp{ik(x — ct)) 
always represent wave propagating at a constant speed c called the phase 
speed. These waves are called progressive waves because they propagate 
without change of form. We know that the role of nonlinearity is to steepen 
a crest while dispersion tends to disperse into waves of different wave- 
lengths. Thus the progressive waves have the role of dynamical equilibrium 
to balance these two effects. 


62 M Rahman 


Normalized Boussinesq equations (41) and (42) for one-dimensional 
waves can be rewritten by defining < u >= дё а8 


дт дф om ; 

8t xm s (n ( x) + = 922 = 0 (47) 
дф дф т? Bo 

ot TES 5 ($2) = 3 8z28t (48) 


Eliminating 7 from (47) and (48) we get a single equation in ф. 


op Fd т 8$ дф 1/89 * 
080 Ər? 3 rt т T *3(8) | (49) 


Rewriting equation (49) in subscript notation for convenience, we ob- 





tain E 
1 
Qu — фы = "g Patt —6 (a T 3%). (50) 
It is this partial differential equation which we want to investigate for 


possible progressive wave solution. We now seek a solution of the following 


form: 
! ф= f) with €=2-ct (51) 
Substituting this form of solution into equation (50), we obtain the 
following ordinary differential equation in f 
(c —1)f derent z)um (52) 
: This equation implies that c = 1 + 0(6,4?). Integrating once with respect 
-to € yields 


(c —1)f = гу" ese (1+) (Р)? + Bi (53) 


But from the first order theory (Airy's theory) we know that n = —¢; = 
cf, and c = O(1) where c = Tae Hence without loss of generality, eqn 
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(53)-can be rewritten as ' 
(c —1)7= TQ t on! + В! | (54) 
It is to be noted that | 
n=- = с, = ccf) р, =cesf') ду". 
Multiplying this equation by т. and integrating with respect to € again 
yields 


! 6 
4-1) = rus tm + Big + В, 


which can be rewritten as 
2 
TT -(- yf ~ 512 — Bin В, (55) 


wher thé integration constants Bı and В» are both of O(6, s^). 
_  Inthe following we shall investigate two important analytical solutions 
of physical interest: (A) solitary waves, (B) cnoidal waves. Part III con- 


cerns with the solitary waves. 


4 Solitary waves 


A solitary wave has a single crest whose amplitude diminishes to zero as 
If] — oo. Thus as || > 00,7 = 1 = тү! = 0 so that В, = 0 and By = 0. 
Equation (55) simply becomes 


п) = 5 ЕЕ - ) т? (56) 


Since the left-hand side of (56) is positive, the right-hand side must 
also be positive for a physical solution. | 

Therefore, we must have c > 1 or in physical variable, С > ./gh, which 
is known as the super critical wave speed. Moreover, у < exl must also ' 


hold good. Hence cot is just the maximum amplitude of the crest which is 


unity because of normalization. Therefore 


—1 
E =1 or с=1+5 (57) 


- 


In dimensional form (57) becomes 


| | 
С = gh (1 + x) or C=[g(h+ A) (58) 


This relation was first found by Rayleigh (1877). Thus substituting 
(57) into (56) yields | 
dn (38) } 
ead Ж ЕР d 59 
d a (Ln (59) 
which can be integrated to give 


G5 (t. = ашата — m (60) 


Thus after reduction it can be written as 
j 
л = sech? (em CON | (61) 


where & is an integration constant and can--be chosen to be zero. In 


dimensional form the surface profile can be written as 


ñ = and E £5 (X — er) (62) 


This is called the solitary wave of Boussinesq (1871, 1872). The solitary 
wave form is shown in Fig. 1. A therefore represents the height of the wave 
and Һ the depth at infinity. 

Solitary-waves can be easily generated in a wave tank by giving any 
kind of impulse. Munk (1949) has done extensive work for applications of 


‘solitary wave theory. Numerical solutions of solitary waves can be found 
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i =0 





Figure 1: Dimensionless free surface profile of a solitary wave 


in the work of Byatt-Smith (1970). Fenton (1972) extended the expansion 
procedure of Benjamin and Lighthill (1954) to the ninth order. In a series of 
remarkable papers, Longuet-Higgins (1973, 1974) and his associates made 
a very thorough investigation of solitary waves. Rahman (1994) recently 
presented a thorough investigation of the interaction of solitary waves. For 
further insight about this topic, interested readers are referred to the work 
of Whitham (1974). 
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On the Pontryagin Maximum Principle for 
Constant — Time Linear Control Systems 
in Locally Convex Linear Topological Spaces 


Hora Krishna Samanta 


Abstract: In this paper some dual characterization of linear systems 
in locally convex linear topological spaces, satisfying the Pontryagin 
maximum principle (briefly, PMP) are obtained. Several examples 
can be found (at least in normed linear spaces) in Rolewicz [7]. 

: Keywords: Locally convex linear topological spaces, constant time 
linear control systems, Pontryagin maximum principle, smooth. 


1 Introduction 


Rolewicz [7] and Singer [9] have discussed constant time linear control 
system in Banach space setting. By a constant time linear control system 
in locally convex linear Housdorff topological vector spaces, we mean & 
system consisting of two real locally convex linear Housdorff topological 
vector spaces X, Y and a continuous linear mapping F : X — Y. A linear 
system is said to satisfy the Pontryagin maximum principle (briefly PMP) 
if for each y(# 0) € F(X), there exists a continuous linear functional 


°AMS Mathematics Subject Classification, Primary : 40118; Secondary : 49J20, 93C15, 
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go € Y" such that 
inf | sup |(2,2^)] : F(x) =y, (2, 2")| < | 
TEX | lex? 
= inf | sup |(z,2^)| : go(F(z)) = golyo), |(т,т')| € if (1) 
TEX TEX’ 
where X*,Y* be the adjoint (conjugate) space to X and Y respectively 
and 0 denotes the null element. 
Our problem is to determine some sufficient condition for a linear sys- 
tem to satisfy PMP in locally convex linear topological spaces over the 


field of real numbers. For this we require the following theorems, and 
definitions. 


Definition 1 We define G^ by 
Gt = {f € X* : f(x) = 0 for all z € G}, (2) 


where G is a linear subspace of a locally convex linear topological 


space X. 
Z 
Theorem 1 The linear system (X LA Y) satisfies PMP iff correspond- 
ing to each ro( 0) € X, there exists go € Y* such that 
sup fiffel: выр Ka) =a} 
f€(ker F)* rcCXCX** 
ICE" (90), 2)] 


Sup | (90), т)| 
ТЄХСХ** 


(3) 


F* denotes the adjoint operator to Ё. 


Proof : Let yo = F(zo). Then, by Hahn-Banach theorem, (Page 21, [5]), 


we get 


| TEX | ex" 


sup — |(zo + KerF,z2’)| 
T'€(X/ Ker F)* 


| t sup |(z,z^)| : F(x) - | 
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= әр (ЛІ: вр 2011). (4) 


f €(kerF)~ 


Again, for апу go(# 0) Є Y*/KerF", the set 


Н = {x € X : |(F'(go), z)| = go(yo) } (5) 


is a hyperplane. Similarly, for any such go, we get, 
inf < sup l(z,z^)|: ggs(F(z)) = 
inf f sap Kx, 2)| : go(F(z)) aoo) 


= sup (20 + KerF"(go),2’)| 
ze(X/KerF'go)* —— 


|(F* (go), zo)! 
вир (Е (ао), 2) ©) 
тЄХСХ*” 


Now if (X 5 Y) satisfies the PMP and zo € X, yo = F'(xo) Æ 0, then, by 
(1); 99(yo) # 0, so F*(go) 5 0; hence by (1), (4), (6) we obtain (3). | 
Conversely if the condition (i) is satisfied and yg € F(X), then, corre- 
sponding to any х0 € X such that yo = F'(zo), from (3), (4), (6) we obtain 
(1). Hence (X E Y) satisfies the PMP, which completes the proof. 
Several equivalent conditions may be obtained in the following way: 


Corollary 1 For a linear system (X EA Y), the following conditions 
are equivalent: 


(i) (X 5 Y) satisfies the PMP; 


(ii) for each zo(# 0) € X with suppcx«|(zo, z')| = 1, there exists 
(F*(go), 2)| = 1 such that 





go € Y* with sup, -ycy 


sup (Iff: mp КА) 


fe(kerF)* 
= |(F* (g0), £0) | (7) 
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(iii) For each zo( 0) € X with 
sup — |(zo + KerF,z)|—1 
=є(Х/КетЕ)' 
there exists go € Y* with 


вир КЕ" (go), 2) = 1, 
rEXCX™ 


such that 
IF" (90), жо)| = 1 (8) 


Corollary 2 If a linear system (X zs Y) satisfies the PMP, it can be 

shown that the image V of F*(Y*) in (X/KerF)*, under the topological 

isomorphism (KerF)~ = (X/KerF)*, is a subspace of characteristic 
І, i.e., 

sup — [(zo + KerF,z’)| 

z'e(X/KerF)* 


= sup|lim- KerF,¥)|: sup ^al-ib — o9) 
wevVv ae 


| zcXcX 
то + KerF € X/KerF, [4] (assuming KerF is a closed subspace of Y ). 


Proof: If the PMP holds, then by F*(Y*) C (KerF)+ and by Theorem 1, 


for each zo € X we have 


sup [Қх0 + KerF,z)| 


a/€(X/KerF)* 

= sup icf: sup M52) 1] 
f €(kerF)* | tEexcx 

= sup filol: mp Ma) =1} (10) 
feF*(Y*) zceXCX" 


Hence, the topological isomorphism f > v of (KerF)* onto (X/KerF)* 
is given by 


\(a + KerF,y)| = |(z, f), 2+ KerF € X/KerF (11) 


we obtain (9), which completes the proof. — \__ 
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Theorem 2 The linear system (X “+ Y) satisfies the PMP, if F(X) 
1з closed in Y. 


Proof: If zo(7 0) € X, then, there exists fy € (KerF)+ such that 


sup (20 + KerF,z^)| = (fo, Zo) | | (12) 
re(X/KerF) вир КЉ, т)| 
reXcX" 


Since F(X) is closed in Y, then (KerF)+ = F*(Y*), by [4]. Bence fy = 
F"(go) for some go € Y*. Using (12) we get (3), by Theorem 1, the linear 
system (X 5 Y ) satisfies the PMP. 


Definition of Polar Set: Let X be a locally convex linear topological 
space and M be a bounded subset of X. Then we shall define the polar set 
M" of the set M C X in X* by (f € X* : вир,„м{[(2, f)|} € 1}, (Page 
136, [10]). | 
Again let E, be the polar set of E C X* in X C X** [10]. Then the 
boundary of E, will be denoted by OE, and is defined as follows: 


OE, = I: € X :supi|(z, f)| = J 
JEE 


Definition of smooth: A locally convex linear topological space X is 
smooth, if at each point x € OE, there is one and only one supporting 
hyperplane of Ey, [2]. 


Theorem 3 Let (X zi Y) be a linear system, with F one to one. 
Then, 


(i) if M" C F*(Y*), then (X 5 Y) satisfies the PMP. 


(it) if X із smooth, the converse also holds. 


Proof: (i) If zo(# v) € X, sup;cx. | (xo, f)| = 1, then, there exists fo € 
(Јо, 2)| = 1, such that |(fo,zo)| = 1, [5]. Hence - 





AX with SUD.excx» 
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fo € М". But, if М" C F*(Y*), then fo = F"(go) for some go € Y". Thus, 
suprexcx~ |(F" (go), 2) = 1, 

I(F* (go), z)| = sup rex: l(zo, f)| = SuPrex/Kerry |(to + KerF, f)| (since F 
is one-to-one). Hence, by Corollary 1, (X Dey) satisfies the PMP. 

(ii) Let X be smooth and that (X — Y) satisfies the PMP, and let 
(fr) = 1. Then by the definition of М", 
there exists то € X with вируех» |(zo, f)| = 1, such that Kro f)| = 1. 
Furthermore, by the PMP and by Corollary 1, for this zo there exists 
go € Y* with sup,exc x» [(F" (go), z)| = 1, such that 


вир |(zo,f) = sup |(ao KerFif)| = |(F" (90), zo)] = 1. 
fex" fe(X/KerF)* 





f € M? with SUDzeXxXcX** 








(since F' is one-to-one). Since X is smooth, we must have f = Е*(9о) Є 
F*(Y*). Hence, since f Є M" is arbitrary with supjexcx~ (f, z)| = 1, we 
obtain by homogeneity, that M* C F*(Y*), which completes the proof. 


Corollary 3 Let (X 5 Y) be a linear system, with X smooth and F 
one-to-one. Then, 
(i) af (X 5 Y) satisfies the PMP, F*(Y*) = X* (in the strong topol- 
ogy); 
(it) af (X zi Y) satisfies the PMP and X is semi-reflexive, F *(Y*) = 
X*, and hence F(X) is closed in Y; 
(iii) If X* is non-separable and F ts compact (X £, Y) does not 
satisfy the PMP. | 


Proof: (i) Clearly M* = X*, by the Bishop — Phelps theorem [1]. Since 

(X 5, Y) satisfies the PMP, (by hypothesis), then, by Theorem 3 (ii), 
М" c F*(Y*). Consequently, X* = M* c F*(Y*), so F*(Y*) = X". 

(ii) Since X is semi-reflexive, so М" = X *. Hence, if (X F, Y) satisfies 

the PMP, then, by Theorem 3 (ii), X* = M* C F*(Y*), so F*(Y*)2xX'. 

- Consequently, F(X) is closed in Y, (Theorem 4, Page 488 [10]). 
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(iii) If F is compact then во is F*, and hence T*(Y*) is separable. But, 
by our assumption, X* is non-separable. Therefore F*(Y*) 5 X*, whence, 
by part (i) above, (X — Y) does not satisfy the PMP, which completes 
the proof. 

We shall give now using Theorem 3(i), an example that the converse of 
Theorem 2 is not valid. | 
Example 1 [2]: Let X = Y - co, where су be the normed linear space 
(co) equipped with the norm ||z|| = sup, |.| where = = £n € со [page 35, 
[10]]. and let 


оо 


‚ Fæ) =F (=) fila)a,, TEE (13) 


i=1 
where, {£n} denotes the unit vector basis of X =- co and {fa} denotes 
the sequence of coordinate functionals on X [so f;(r,) = 6;;]. Then F is 
one-to-one and 


oc 


Р") = Y Güg(z)fs 9 €Y* = х" (14) 


1=1 


$ 
on the other hand, it is well known (вее Remark (+) to Theorem 5, [6]) that 
М" = lin{ fn}, 


the linear span of {f,}. Thus by (14), X* = M* C F*(Y*), whence, 
by Theorem 3(), (X £, ү) satisfies the PMP. However, by (13), F(X) = 
X =Y, F(X) z X,so F(X) is not closed, which shows that the converse of 
Theorem 2 is not valid. Note also that , in this example, Ё' is compact, but 
X" is separable and X is not smooth and not EE ув [see Corollary 
3 (ii) and (iii) above}. 
‚Now F is one-to-one in Theorem 3, Corollary 3, and Example 1 and the 
fact that F(X) = Y in Example 1 are not accidental, since the following 
" theorem powe, that the general'case can be reduced to the case when F 


has both of these properties. 
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Theorem 4 For а linear system (X 5 Y), the following statements 
are equivalent: 


(i) (X 5 Y) satisfies the PMP 
(ii) (X/KerF = Y) satisfies the PMP where F is thé (one-to-one) 
mapping induced by Е, that 15 


F(x + KerF) = Е(х), = + KerF є X/KerF; (15) 


(ii) (X 5 F(X) satisfies the PMP. 


Proof: By ‘Theorem 1, we have (it) iff, for each zo + KerF є X/KerF, 
there exists ду € Y* such that 
sup Цао + Ker): sup Iz) 1] 
Ve(KerF)" rexXcx"- 
= |( + KerF,P*(go)): sup |(F*(90),2)| (16) 
|"r€XCX" 


< But, since F is one-to-one [indeed, F(x + KerF) = 0 implies, by (15), 
x € KerF, whence x + KerF = Ker F, the zero element of X/KerF], we 
have (KerP)- = (X/KerF)*; hente; by the canonical linear isometry (11) 
of (KerF)+ onto (X/KerF)*, the left hand sides of (3) and (17) coincide. 

Furthermore, by (15), we have 


|" (zo + Ker F), 90)| 
|(zo, F* (go))| (17) 


so F*(go) is the image of F*(go) by the canonical linear isometry (11) of 
(KerF)+ onto (X/KerF)*, and therefore 


sup |(F*(go),z)|- вир |(F*(gq),2)| 
rcXcx* rcexcx" 


(х0 + KerF, F*(go))| 
= |(Р(20), 90)| 


This together with (17), shows that the right-hand side of (3) and (16) 
. coincide, too. Thus (i) Ф (ii). 
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у 
Finally, the equivalence (i) <= (iii) follows from the definition (1) of the | 
PMP and the observations that go | F(x)€ (F(X))* for each gy € Y* and 
that, conversely, each үр € (F(X))* can be extended to a œ € Y*. This 
completes the proof of Theorem 4. 
Combining Theorem 4 and Corollary 3 (i) and (it), we obtain the fol- 


lowing generalization of Corollary 3(2) and (iż). 
Corollary 4 Let (X 5 Y) be а linear system, with X smooth. Then 
` (i) if (X SY) satisfies the PMP, F*(Y*) = (KerF)!; 
(ii) if (X ds Y) satisfies the PMP and X is semtreflezive, F*(Y*) = 
(КетЕ) and F(X) is closed in Y. 
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Flow of a Bingham Fluid 
by means of Peristaltic Waves 


VP Srivastava 


Abstract: The problem of peristaltic induced flow of a Bingham 
fluid under long wavelength approximation is investigated. It is found 
that the pressure drop increases with the flow rate and there exists a 
linear relation between them. For any given set of other parameters, 
the magnitude of the pressure drop is higher for Bingham fluid than 
Newtonian fluid. Pressure drop decreases with increasing amplitude 
ratio but increases with the yield stress. The friction force possesses 
the character similar to the pressure drop. 


1 Introduction 


Fluids transport by means of peristaltic waves has' been the subject of - 
engineering and scientific research since the first invstigation by Latham [1]. 
It is a form of fluid transport that occurs when a progressive wave of area 
contraction or expansion propagates along the length of a distensible tube 7 
containing liquid or mixture. Besides its applications to physiological flows 
[2], the study of the mechanism of peristaltic pumping is of significant use to 
engineering works. To understand peristaltic action in different Situations, 


several theoretical and experimental studies have been conducted in the 
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literature. An account of the most of the investigations reported so far 
may be had from Srivastava and Coworkers [2,3]. 

Excepting a few, most of the studies in the literature have been carried 
out for Newtonian fluids. It is known that a large number of fluids of 
practical importance including the physiological fluids behave like non- 
Newtonian fluids. The studies [2-7] have considered this aspect of the 
problem since the initial investigations by Raju and Devanathan [8,9]. The 
research reported here considers the peristaltic transport of a Bingham 


fluid. a special case of a non-Newtonian fluid which inhibits yield stress. 





Fig 1 Geometrical representation of peristaltic waves In 
a circular tube.. 


2 Formulation and Analysis 


Consider the axisymmetric flow in a circular tube with a sinusoidal wave 
travelling down its wall. The geometry of the wall surface is described as 
(Fig. 1) 
‚ 2л 
H (x, t) = а + bsin (x — et), (1) 


where a is the radius of the tube, b is the amplitude of the peristaltic wave, 


- 
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А is the wavelength, c is the wave propagation velocity, t is the time and x 
is the axial coordinate. 

Assuming that the circular tube is filled with à homogeneous and in- 
compressible Bingham fluid, the appropriate equation in the wave frame 
of reference (moving with speed c) under long wavelength approximation 
and neglecting the inertia terms [10], is expressed as 


19 dp 





-a M irr] — — 2 
(тты) = —5 (2) 

where TT,r, for a Bingham fluid. is given by 

=a МЕЕ Ты Tos 

ти m Ts ct T (3) 
Ou 0 2 (4) 

uu = ; Tp , 

Or £e 


where р is the pressure, p is the velocity of the fluid, тр is the yield stress, r 
is the radial coordinate and и is the axial velocity of the fluid. The relation 
(4) corresponds to the vanishing of the velocity gradient in the region in 
which Tys < ту and implies a plug flow. ` 


Ап introduction of the following dimensionless variables 


ct 
r= Е г! = t из 
и = = TQ = p р = LP (5) 
TH equations (2)- (4), after dropping primes, yields 
1с. [ E- +n) = -Ф, Тек = To, (6) 
= = 0; Tre < To, (7) 


where Tyr and To are dimensionless shearing and yield stresses respectively. 


The non-dimensional boundary conditions are 


u = —l at r = H/a-h-1-désin2mz, (8) 
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7,,18 finite аёт = 0 to ap, (9) 


where ф = b/a < 1 and ар is the radius of plug flow region. 
The straight forward integration of equation (6) under the boundary 


conditions (8) and (9), yields the expression for velocity, u as 


«=-ї+гы (525 tn) - (22 + ») |. (10) 


The constant plug flow velocity u = up is obtained from equation (10) 
evaluated at r = ар, a, being the non-dimensional radius of plug flow 
region and is given by [11] 
dp ay 
———— 11 
n=- 29 п) 
which: may be conveniently calculated using the expression for the flow 
rate. 
Considering the radius of plug flow area to be small as compared to the 
rion-plug flow area [3,13,14], one obtains the expression for the instanta- 
neous volumetric flow rate Q (= Q'/ na*c, Q' being the flux in the moving 


system) as 


& 
| 
N2 
om 
> 
з 
e 
ti. 
ш) 


1 
~ 
bo 
+ 
су“ 
Aa 
- 
mA 
гы, 
ке 
| 
SIG 
eee” 
t 
=ч 
Y 
-+ 
Co] bo 
SIG 
a | = 
PETN 
SIE 
ee” 
bo 


4 
a2 (3 |) (12) 
Oo kT i 
where m, = (—h/2)(dp/dx) is the shear stress at the tube wall. Under the 
condition [13,14], 70/74 « 1, equation (12) yields 


x 2 1 ту 
xz (Qe). (13) 
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Now following the report by Shapiro et al. [10], one derives the mean 


volume flow rate, Q over a period as 


g-Q«i« f. (14) 


An use of equations (13) and (14) into equation (11) determines the 
radius of plug flow region in its non-dimensional form as 
E 3mh’ 
P 4B(Q-1- g? + №2) + nh] 
which reveals that the plug flow radius depends on yield stress, the tube 


(15) 


geometry and the flow rate. a fact which has been previously observed in 
[2,12,15] etc. 

Since the pressure drop, Ap = p(0) — р(А), across one wave length is 
same whether measured in moving or stationary coordinate system, it is 


therefore calculated using equation (14) into (13) as 


1 dp 
d 

8 1 = 3 
пут Е - фт + QU + 39) 


9 gn 
+ 5-04-10). (16) 


Ap 


The dimensionless friction force per wavelength, F (— F ' IrAcp, F' being 
the friction force and is same whether measured in moving or stationary 


system) at the wall is given as 


e 
=- зи 0-1-5 +(1+2) ace. (17) 


The pressure drop for zero time-mean flow, (Ap)g= and the time-mean 


F 


| 


flow for zero pressure diop, (Q)Ap-o which are of particular interest, from 


A 
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equation (16) are obtained as 


(Ар) = : La — dn + (9 —16 
Pg= = (gym 3! -$)Tn-tr-(ó—-160|, (18) 


(Quap-o = gag [^ (-$) «se - vn), m 


In absence of the plug flow (i.e., ту = 0), equations (16)-(19) yield the 


corresponding results for a Newtonian fluid as 


Ap = йй ба E: 59) + Pe Е z ' (20) 


8 = ф? 
t сег ро __ A248/2 

25и |016 aon), om 

29^ (d* — 16) 
(Ap) 5 1 ey” 
$^(16 — 9?) 
2(2 + 342) ` 
Equations (20)-(23) are the same results as derived from the analysis of 
Shapiro et al. [10], Jaffrin and Shapiro [16] and Barton and Raynor [17]. 


(22) 


(Q) apo (23) 


3 Numerical Results and Discussion 


In order to discuss the results of the analysis quantitatively, computer 


.codes were developed for the numerical evaluations of the analytical results 


obtained for dimensionless pressure drop and the friction force (equations 
(14) and (15)). Some of the important results are displayed graphically in 
Figs. 2-6. 

The average pressure drop, Ap versus the time-mean flow rate, Q is 
plotted for different values of the parameters ф and ту, which shows а 
linear relation between them. As expected, an increase in the flow rate 


increases the pressure drop and thus minimum flow rate is achieved at zero 
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Figure 3 Variation of pressure drop, Ap with Ø for different Q and Co 
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pressure drop and the minimum pressure drop occurs at zero flow rate. The 
magnitude of the pressure drop under a given set of parameters is greater 
for a Bingham fluid than for a Newtonian fluid. The variation of pressure 
drop Ap with amplitude ratio ф is shown in Fig. 3. The pressure drop 
decreases indefinitely with increasing amplitude ratio for any given set of 
other parameters. One notices from.Fig. 4 that pressure drop, Ap increases 
with the yield stress, 7 and for any given set of other parameters, its 
minimum value is achieved at zero yield stress (Newtonian fluid). A study 
of Figs. 5 and 6 reveals that the friction force, F possesses a character 
similar to the pressure drop, Ap (an opposite character to the pressure rise 
—Ap). 
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